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Abstract: We consider the functional Schrodinger equation for a self interacting scalar 

field in an expanding geometry. By performing a time dependent scale transformation on 

. the argument of the field we derive a functional Schrodinger equation whose hamiltonian is 
O . 

■"sj- , time independent but involves a time-odd term associated to a constraint on the expansion 



current. We study the mean field approximation to this equation and generalize in this case, 



^ for interacting fields, the solutions worked out by Bunch and Davies for free fields. 
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Several authors have pointed out that inflationary models based on effective potentials ignore 
the fact that the very concept of effective potential implies a static approximation which may 
be invalidated by dynamical effects 0. The purpose of this present letter is to investigate this 
question using the functional Schrodinger equation for a scalar field in an expanding geometry: 

= Ij^—J^g^ + Wt )(V*>(x)r + a''M( m ^( X ) + A^ W )]» } . (1) 

In this equation a(t) = exp(xt) specifies the metric ds 2 = dt 2 — a 2 (t)<ix 2 and \ is the expansion 
parameter (Hubble's constant). The constants m and A are respectively the bare mass of the 
field and the bare coupling constant. The evolution equation ([[]) was considered by Guth and 
Pi || and by Eboli, Jackiw and So- Young Pi 0. In particular in reference || a numerical 
solution of this equation was attempted in the mean field approximation. However instabilities 
were found in the implementation of the method. Our aim in the present study is to perform 
a particular scale transformation on the wave functional which brings equation ([I]) into the 
form of a static functional Schrodinger equation. The method, although more elaborate, is 
reminiscent of the transformation used in nuclear physics to describe rotating nuclei |Q. 
Let us look for solutions of equation flU) in the form 

$[^x),t] = vI/[^(x),t]^L=, (2) 

'N(a) 



where 

£(x) = <p(a(t)x), (3) 

and where N(a) is a normalisation factor. By using the definition of the functional derivative 
we first find 



5$ 



5^ r 5$ 

— ^(ax)rfx = / —— 

which allows one to express the functional derivative of $ in terms of the functional derivative 
of * as 

5$ 1 5^ 



<5y?(x) a 3 5£(x/a:) 
By applying this formula twice we find that the kinetic energy of the field is 



(5) 



8 2 § 1 r 5 2 V 



The transformation law for the mass term is easily derived from the relation 

J dx^ 2 (x) = a 3 J dy<p 2 (ay) = a 3 J dy £ 2 (y), (7) 

while for the gradient term we have 

/ rfx (V^(x)) 2 = ^ J dy (V y ^(av)) 2 = a J dy (V£(y)) 2 . (8) 

We still have to calculate the time derivative of the wave funtional $. Between t and t + dt the 
change <5$ in $ is obtained by combining the time change of \l/ and the change in a(t) with the 
result 



Since the change in £ is 



5ct 5cx 

<Jf (x) = tp(ax + x5a) - <p{ax) = f (x H x) - f (x) = — x.V£ , (10) 

a a 

we obtain 

a r N 
d t $ = d t m + i-J dx7r(x)x.V£^- — * . (11) 

The time derivative of N(a) is easily calculated in a similar way. The final result is 

d t $ = d t y + -- / dx{7r(x)x.V£ + x.V£vr(x)}^ . (12) 
2 a J 

By choosing a(t) = l/a(t) = exp(— %t), we thus see that $ is a solution of the time-dependent 
Schrodinger equation ([!]) provided \1/ satisfies 

zd t m = H x m, as) 

where H x is the time- independent hamiltonian 



<5£(x)<S£(x) 

The last term in the previous equation is static. However it corresponds to a constraint on the 
expectation value of the expansion current. Indeed in ordinary quantum mechanics of a single 
particle this operator would be similar to the operator x.p + p.x. 



The presence of this term illustrates the relevance of the discussion given in reference KJ 
about the importance of dynamical effects. Indeed, as will be shown below, this term affects 
the gap equation and may thus drive the system out of the symmetric or asymmetric phase. 

In order to obtain a mean field approximation to the functional Schrodinger equation (13) 
we assume that at each time t the wave functional describing the system is a gaussian functional 
of the form: 



ATexp{>/^(x,t) («x)-fix,t) 
xexpj-jdxdy (£(x)-e(x,0) (^(x, y) + zE(x, y, *)) (£(y) - £(y, *)) } 



+ ! 



(15) 

where M is a normalization factor. Our trial functional contains four variational parameters: 
the center of the gaussian £(x, t), its conjugate momentum 7f(x, t), and the real and imaginary 
parts of its kernel: G(x, y, t) and E(x, y,t). 

The expectation value of the hamiltonian H x evaluated with the gaussian trial state \1> is 

< m x \* >= i ^{itt 2 + kvo 2 + \mie + ±\e 

"jG-^x, x, t) + 4 (GEG) (x, x, t) + (-A + m 2 + \? + ±AG(x, x, £)) G(x, x, t) 

"7f(x,t)x.V^(x,t) -£(x,t) V.X7f(x,t)] -xTrace{S(x.V)G-G(V.x)S}} . 

(16) 

By varying the action || < \l/|i<9i — H x \§> > with respect to the variational parameters, we 
obtain the following dynamical equations : 

d t G= 2(GS + EG) -x{(x.V)G-G(V.x)}, 

d t E = |G- 2 - 2£ 2 - §{-A + ml + ±Af + ±A G(x, x)} 

- X {(Ex.V)-E(V.x)} (17) 
dtl= -7f + x(x.V0, 

a t 7f = x(x.Vvr) + {-A + m 2 + ±A£ 2 + \\ G(x, x)}£. 

To derive these expressions one has to remember that the cyclic invariance of the trace does 
not hold for pairs of operators such as x and p. 
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The previous equations admit translationally invariant solutions of the form 

e(x,*)= m, 

7f (x, t) = 7f (t) , 

(18) 

G ( x ' y) = (i^F ^ rfk G ^ ex p( ik ( x - y)) > 
s ( x > y) = (2^f ^ dk E ( k "> ex p(^ k ( x - y)) ' 

provided that £, n, G(k) and satisfy : 
d t G{k) = 4G(k)Z(k) - X {(k.V fc )G(A;) + 3G(k)} , 

^S(fc) = lG-\k) - 2^(k) + |{k 2 + ml + §A£ 2 + | ' rfkG ( fc )i " X{3S(^) - k. V fc S(A;)} 
ftTf = (rug + \Xl 2 + f / ^k G(*)) £ • 

(19) 

Static solutions in the expanding frame correspond, after elimination of tt and £ to : 



m » + ^ + I(2^/ AGW }^=° (20) 



and 

^ ^ + g " *I g - 5 g ~ 2 * 2 (I g ) 2 - w + 2 (5 - - 2 ) ^° • (21) 

where v satisfies the self-consistent equation : 

„ Y ^_ mS _^_y_/ dkG(k ). (22) 

Equation fl2T|) is the differential equation for the square root of the modulus of the Hankel 
Function, and the most general solution has the form: 

G(k) = I- Cl H^{k/ X ) + c 2 H^\k/ X ) 2 (23) 

where are the first (i=l) and second (i=2) Hankel function of order u, and the coeficients 
Ci and C2 obey the relation: |ci| 2 — | C2 1 2 = 1, and are determined by the boundary conditions. 
If we impose that the asymptotic behavior be the same as in flat space-time we get C\ = 1 and 
c 2 = H. 



Equation (|22|) will be refered to as the gap equation since it generalizes the equation obtained 



for the mass gap of a self interacting scalar field in mean field theory in Minkowski space. This 



limiting case is obtained for x = 0. In the limit \ = 0, v is large and imaginary and behaves 
as im/x where m is the mass gap at x = 0. Using the asymptotic formula 

* ^L- , (24) 

one obtains = 1/2^/k 2 + m 2 as expected. 

Inserting the gap equation into the equation for the background field £ we find that static 
solutions are determined by the two conditions : 



"Y - ~ A X 2 + ml + ^ + [ dk F« 2 = (25) 



and 



(y 2 x 2 + \t 2 - \x 2 ) t = o. (26) 



Before discussing the properties of these two equations it is worthwhile examining the structure 
of the associated solutions when they are rewritten in the laboratory frame. Denoting by G 
the kernel of the gaussian in this frame we have G(k) = e~ 3xt G(ke~ xt ), i. e. 

G(k,t) = e-^^- Hl(ke~ xt ) 2 . (27) 

For a free scalar field, v is given by 

v 2 X 2 = \x 2 -™l (28) 

and the expression (B^) is identical to the solution obtained by Bunch and Davies 0. In our 
formalism this solution appears as a static solution in the expanding background geometry. 
However, our solution ( p7|) holds also for the case of a self interacting scalar field. Different 
classes of solutions with broken or unbroken symmetry can be found. 

As for the case % = the gap equation contains divergent quantities and must be renor- 
malized. In order to investigate this question let us introduce a momentum cutoff A. At large 
momentum, using the asymptotic form of Hankel functions, one finds that G(k) behaves as 
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so that the quantity G contains a quadratic and a logarithmic divergence when A goes to infinity 
^)^{A-^ l0g (^_)}. (30) 

The key to the renormalization of our variational calculation is the existence, for a given value 
of A, of a critical value r c of the ratio r = m^/A 2 beyond which there is no asymmetric solution 
of the gap equation 0. For a large value of the momentum cutoff, the gap equation can indeed 
be written as 



2 ' , X 2 2 , \ F 2 , Ar 2 2 X 2 \ , fX — 4l/ 2 X 2 . 



where e = ±1 specifies wether r is above or below r c and where we have introduced the 
following quantities 

2 + log ($r) 



and 



V 2 2A fl A 2 X \ 

e A^ = ^r (ro + To^- a*-)' (33) 



i.e. 

4 = ! r (r »- r J- (34) 

Keeping the ultraviolet cutoff A fixed, the approach to the critical point depends on the value 
of the expansion parameter \. 

Below r c (i.e. e = — 1) there is a solution with broken symmetry while for ro > r c there is 
only a symmetric solution 0. For a given value of A, approaching the critical point from the 
left the solution of the gap equation is 

2 2 , X 2 4 2 _16 ? £ 

-^X+Y = - A e * . (35) 

We thus see that when A goes to infinity this solution remains finite provided that the bare 
coupling constant A evolves as 

' 1 logf^-V (36) 



A 16vr 2 \M 2 R j' 



where Mr is an arbitrary scale. 
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In conclusion, we have been able to find self similar solutions of the mean-field evolution 
equation for the fluctuations of a self interacting scalar field in an expanding geometry. This 
has been possible in the particular case of a de Sitter geometry for which the expansion rate 
a/a is a constant. In the case of a radiation dominated univers for which a(t) = y/t, this type of 
solutions can not longer be found. However, in this case one has a/a — l/2t and the formalism 
we have described can be used to make predictions about the long time behaviour. 
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